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Abstract 

We give an overview over models in which cosmic rays above ~ 1 EeV (= 10 18 eV) 
are produced by the decay of supermassive "X" particles released from topological 
defects possibly created in cosmological phase transitions. We note that, for an 
interesting particle physics parameter range, these models are still consistent with 
current data, and discuss signatures for the topological defect mechanism which 
can be tested by the next generation experiments. 



1 Introduction 

The highest energy cosmic ray (HECR) events observed above lOOEevSH are 
difficult to explain within conventional models involving first order Fermi accel- 
eration of charged particles at astrophysical shocks. On the one hand, even the 
most powerful astrophysical objects such as radio galaxies and active galactk: 
nuclei (AGN) are barely able to accelerate charged particles to such energies. u 
On the other hand, above ~ 70 EeV the range of nucleons is limited by photo- 
pion production on the cosmic microwave background (CMB) to about 30 Mpc 
(which is known as the GZK-effectlj), whereas heavy nuclei are photodisinte- 
grated on an even shorter distance scale. □ In addition, for commonly assumed 
values of the parameters characterizing the galactic and extragalactic magnetic 
fields, protons above f 00 EeV are deflected by only a few degrees over these 
distances a and obvious powerful sources have not been found in the arrival 
direction of the observed HECR events. 13 

Within "top-down" (TD) scenarios, in contrast, predominantly 7-rays and 
neutrinos are initially produced at ultra-high energies (UHEs) by quantum me- 
chanical decay of supermassive elementary X particles related to some grand 
unified theory (GUT). Such X particles could be released from topological 
defect relics of phase transitions which might have been caused by sponta- 
neous breaking of GUT symmetries in the early universe.! Topological defects 
from phase transitions in the early universe such as cosmic strings, monopoles, 
and domain walls are topologically stable, but nevertheless can release part 
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of their energy via collapse or annihilation in the form of X particles. The X 
particles can be gauge bosons, Higgs bosons, superheavy fermions, etc. depend- 
ing on the specific GUT, their mass mx being comparable to the symmetry 
breaking scale. They subsequently typically decay into a lepton and a quark 
which roughly share the initial energy equally. The quark then hadronizes into 
nucleons (Ns) and pions, the latter ones in turn decaying into 7-rays, elec- 
trons, and neutrinos. Given the X particle injection rate, dnx/dt, the effective 
injection spectrum of particle species a (a = j,N,e^,v) can be written as 
(dnx/dt)(2/mx)(dN a /dx), where x = 2E/mx, and dN a /dx is the relevant 
effective fragmentation function. We take the primary lepton to be an elec- 
tron injected with energy mx/2. The total hadronic fragmentation function 
dNh/dx can be taken from solutions of the QCD evolution equations in modi- 
fied leading logarithmic approximation which provide good fits to accelerator 
data at LEP energies. u Fig. [j] shows this function-for mx — 2 x 10 16 eV in com- 
parison to approximations used in earlier work. 113b Motivated by LEP data at 
lower energies, we assume that about 3% of the total hadronic content consists 
of nucleons and the rest is produced as pions and distributed equally among 
the three charge states. The standard pion decay spectra then determine the 
injection spectra of 7-rays, electrons, and neutrinos. The X particle injection 
rate is assumed to be spatially uniform and in the matter-dominated era can 
be parametrized El as dnx/dt oc t~ 4+p , where p depends on the specific defect 
scenario.. The case p = 1 is representative for a network oLprdinary cosmic 
strings til and annihilation of monopole-antimonopole pairs. Ii3 

Sincq-the absolute flux level predicted by TD models is very model de- 
pendent, tia we will allow an arbitrary flux normalization noting that certain 
TD scenarios such as annihilation of magnetic monopole-antimonopole pairs£3 
yield HECR fluxes consistent with observations without violating bounds on 
monopole abundances. Such models are attractive in explaining HECRs be- 
cause they predict injection spectra which are considerably harder than shock 
acceleration spectra and, unlike the GZK effect for nucleons, there is no thresh- 
old effect in the attenuation of UHE 7-rays whose mean free path in the cosmic 
low .energy photon background is probably larger than that for nucleons (see, 
e.g.M). 

2 Signatures for and Constraints on TD Models 

Based on the general features of the type of TD scenarios discussed in the pre- 
vious section, there are two distinctive signatures for them at energieSpabove 
100 EeV: First, the observed primaries should be predominantly 7-rays. Sec- 
ond, there could be a spectral feature in the form of a "gap" .113 Increase of the 
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Figure 1: Approximations to the hadronic fragmentation function dN^/dE for rax = 2 X 
10 16 GeV. The solid line indicates tie solution to the QCD evolution aquations in modified 
leading logarithmic approximation, □ the dashed line is the one used inrfj p*ifl-*hc dotted line 
represents the approximation dN h /dx = (15/16)x~ 1,5 (1 — x) 2 used infcj' Ha and in Fig. 0. 



current total exposure at these energies by factors of a few could distinguish 
between acceleration type sources and TD mechanisms at the 99% confidence 
level. 113 This should easily be possible with next generation experiments un- 
der construction such as the High Resolution Fly's Eye E_ZJ or in the proposal 
stagelia such as the Pierre Auger project. 

Recently, there has been considerable discussion in the literature whether 
the 7-ray, nucleon, and neutrino fluxes predicted by TD sce narios a re consis- 
tent with observational data and constraints at any energy. Ejc3'EJ To thor- 
oughly investigate this, we have performed extensive numerical simulations 
for the propagation of extragalactic nucleons, 7-rays, and electrons with en- 
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ergies between_L0 8 cV and 10 25 eV through the universal low energy photon 
background, oE3 which includes the radio background, the CMB, and the in- 
frared/optical (IR/O) background. All relevant interactions have been taken 
into account, including synchrotron loss in the EGMF of the electronic com- 
ponent of the electromagnetic cascades which result from UHE 7-ray injection 
into the universal radiation background. 




E (eV) 

Figure 2: Predictions for the average differential fluxes of 7-rays, nucleons and muon and elec- 
tron neutrinos by a typical TD scenario as described in the text. Vv4s used the hadronic frag- 
mentation function in modified leading logarithmic approximationta for mx = 2 X 10 25 eV. 
The average EGMFrstrength was assumed to be 10 — 12 G. Also shown are the combined data 
from the Fly's Eye El and the AG AS A □ experiments above lOEeV (dots with error bars), 
piecewise power law fits to the observed charged KR flux (thick solid line) and experimental 
upper limits on the 7-ray flux below 10 GeV fromEJ (dotted linereji left margin). The arrows 
indicate limits on the 7-ray flux from. £3 

Fig. ^| shows the results for the 7-ray and nucleon fluxes from a typical TD 
scenario, assuming an EGMF of average strength 10~ 12 G, along with current 
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observational constraints on the 7-ray flux. The spectrum was normalized 
in the best possible way to allow for an explanation of the observed HECR 
events, assuming their consistency with_a nucleon or 7-ray primary (although 
a 7-ray primary is somewhat disfavoredES) . The flux below a few tens of EeV is 
presumabhz_caused by conventional acceleration. The likelihood significance of 
the fit (seePJ for details) in Fig. | is ~ 50% for the energy range above TOO EeV. 
While the shapes of our spectra are similar to the ones obtained in, Eil this is 
in contrast to their procedure of normalizing to the observed differential flux 
at 300 EeV which overproduces the integral flux at higher energies. 

Since the comparatively large amount of energy injected at high redshifisis 
recycled to lower energy 7-rays, TD models are significantly constrained P^H 20 ! 
by current limits on the diffuse 7-ray background in the 100 MeV — fOGeV 
region. Ej Note that the IR/O background strongly depletes the 7-ray flux in 
the range fO 11 — 10 14 eV, recycling it to energies below fOGeV (see Fig. ||). 
Constraints from limits on CMB distortions and light element abundances 
from 4 He-photodisintegration are comparable to the bound from the directly 
observed 7-rays. c3 The scenario in Fig. || obeys all current constraints within 
the normalization ambiguities and is therefore quite viable. 

Whereas the UHE nucleon and 7-ray fluxes are independent of cosmo- 
logical evolution, the 7-ray flux below ~ f 11 eV and the neutrino flux are 
proportional to the total energy injection which, for all other parameters held 
fixed, increases with decreasing p. Ej For mx — 2 x 10 25 eV, scenarios with 
p < 1 are therefore ruled out (see Fig. ||), whereas constant comoving injection 
rates (p — 2) are well within the limits. Since the EM flux above ~ f 22 eV is 
efficiently recycled to lower energies, this constraint turns out to be basically 
independent of nix in case of a low EGMF,E3 in contrast to earlier analytical 
estimates based on the continuous energy loss approximation which underesti- 
mates the flux around f 00 EeV.oEj The constraints from the flux limits below 
lOGeV become somewhat tighter for an EGMF of strength > f0 -11 G. 

The predicted neutrino fluxes @ are also consistent with bounds from the 
Frejus experiment. c3 At these flux levels, neutrinos are unlikely candidates 
for the observed-IJECR events due to their small interaction probability in 
the atmosphere. E3 A future detection of an appreciable neutrino flux above 
~ 1 EeV, for example, by a km 3 scale neutrino observatory cB could establish 
an experimental lower limit on the ratio of energy injected as neutrinos versus 
electromagnetically interacting particles and thus probe GUT scale physics. c3 
In the scenario shown in Figs. [I] and this ratio is about 0.3. 

Our simulations show that an isotropic 7-ray to total CR (7/CR) flux ratio 
at ~ 10 EeV as high as 10% can be attained. However, this is only possible 
if a TD mechanism of the type discussed above is responsible for most of the 
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HECR above ~ 100 EeV, and if the EGMF is weaker than ~ lO" 11 G on scales 
of a few to tens of Mpc. In case of acceleration sources of the HECRs, such 
high 7/CR flux ratios at ~ 10 EeV can only be attained in the direction of 
powerful nearby acceleration sources. In contrast to the TD models, in this 
case the 7-rays would be produced as secondaries of nucleons interacting with 
the CMB. The spectral shape of the 7-ray flux also depends on the EGMF 
which determines where the energy loss of 7-rays is dominated by synchrotron 
loss rather than inverse Compton scattering on the CMB and_thus could be 
used to "measure" the EGMF in the range ~ 10~ 10 - 10~ 9 G.E3 

Information on the EGMF structure could also be obtained by observing 
the energy and and arrival time distribution of nucleons from sources which 
release UHE^psmic rays on a time scale short compared to ~ 1 yr, i.e. in 
a burst. E2lEilt3. This is because the average time delay caused by deflection 
in the EGMF is ~ 0.9(£:/100EeV)" 2 (r/10Mpc) 2 ( J B/10- 11 G) 2 (7 C /1 Mpc) yr, 
where E is the observed nucleon energy, r is the source distance, and B and 
l c is the typical strength and the coherence length of the EGMF, respectively- 
Some TD scenarios such as the ones involving collapsing cosmic string loops liil 
would imply the existence of such bursting sources. 



3 Conclusions 



Some TD type scenarios of HECR origin are still unconstrained by current 
data and bounds on 7-ray and UHE CR fluxes. For example, in case of an 
EGMF < 10 -9 G, spatially uniform annihilation of magnetic monopoles and 
antimonopoles is still a viable model for GUT scales up to 10 16 GeV. In such 
scenarios, the flux would be dominated by 7-rays above ~ 100 EeV and the 
possibility of a gap in the spectrum arises. A solid angle averaged 7/CR flux 
ratio at the 10% level at ~ 10 EeV is a signature of a non-acceleration origin 
of HECRs hinting to the presence of a TD mechanism. At the same time it 
would put an independent new upper limit of ~ 10~ n G on the poorly known 
EGMF on scales of a few to tens of Mpc. 

The fact that some TD mechanisms would imply bursting sources could 
provide another "probe" of the EGMF. 
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